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We study driven Id lattice gas models with two types of particles and nearest neighbor hopping. 
We find the most general case when there is a shock solution with a product measure which has a 
density-profile of a step function for both densities. The position of the shock performs a biased 
random walk. We calculate the microscopic hopping rates of the shock. We also construct the 
hydrodynamic limit of the model and solve the resulting hyperbolic system of conservation laws. In 
case of open boundaries the selected steady state is given in terms of the boundary densities. 
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I. INTRODUCTION 

O ' On the macroscopic level driven diffusive systems are often described by hydrodynamic equations for some relevant 
parameters (usually the particle densities). These partial differential equations are generally nonlinear and can develop 
shocks, i.e., discontinuities in the space dependence of the densities. These shocks behave as collective excitations in 
systems with one conservation law: they can be characterized by only one parameter, namely their position, and the 
propagation can be described by single particle dynamics. 

Recently much attention was payed to the investigation of the microscopic structure and the microscopic dynamics 
of such shocks PI 0, 01 • Mostly the well-known asymmetric simple exclusion process (ASEP) was studied: it 
was pointed out in U (and for infinite systems in Q) that for special tuning of densities and microscopic hopping 
rates there exists a travelling shock with a step-like density profile even on the microscopic scale, which behaves like 
a one-particle excitation. However, little is known about the microscopic structure of shocks in systems with two 
conservation laws |E0)H>@) which recently have become a focus of attention (for a review see [Io|). 

In the present work we show that a shock measure with single-particle dynamics can describe also systems with 
two conserved densities, introduced in section [H] The position of the shock performs a biased random walk just like 
in the ASEP (section ITTT|l . In section llVl we also study the hydrodynamic limit of the model (under Eulerian scaling) 
which shows a larger class of stable shock solutions. The hierarchical structure of the hydrodynamic equations for 
this system allows us to deduct the steady state selection in an open system connected to particle reservoirs at its 
boundaries. 

1 1 II. TWO-SPECIES EXCLUSION PROCESS 

O 

The model is defined on an open lattice with L sites and two types of particles (A and B). Each lattice site can 
i either be occupied by a particle, which can be either A or B, or can be vacant (0). Having two independent conserved 
quantity requires that only hopping processes are allowed. For simplicity we allow just for nearest neighbour hoppings. 
£h These are the following: 

■ A0 -> 0A 0A^ A0 

Ch ! B0^ 0B 0B^B0 

q ! BA — > AB AB -> BA (1) 

A state of the model is defined through a probability measure P^ on the set of all configurations rj = (771, r)2, • ■ ■ ,T]l), 
r]k € {A,B,0}. For our purposes it is convenient to use the Hamiltonian formalism where one assigns a basis 
vector \r]) of the vector space (C 3 )® L to each configuration and the probability vector is defined by \P) = J2 V Pvl 7 ))- It 
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is normalized such that (s\P) = 1 where (s\ = (^l- The time dependence is now described by the master equation 

\P(t)) = -H\P(t)) (2) 



J7) \ 

d 

df' 



through a "Hamiltonian" H which has matrix elements H^^i the hopping rates between configurations 77, rf . Since 
we have only nearest neighbor exchange processes the Hamiltonian in J5J can be written as 



L-l 



H = hi + fefc,fc+i + h L , (3) 



fc=i 
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where h^^+i acts nontrivially only on sites k and k + 1 (corresponding to hopping) while h\,hL generates boundary 
processes specified below. Let If., 1 ., 2 (0i,0 2 ,0i,02 e {A-B> }) be an operator on C 3 ® C 3 which acts on the basis 
vectors 1 771,772) as 



W°$\ m ,r, 2 ) = 6 m o 1 8 V2e2 \d' 1) 6' 2 ). 



Then hk,k+i can be written as 



_ n®(fc-i) 



E 



^1^2 V ^2 



(L-fc-1) 



(4) 



(5) 



The diagonal term in (J5J stands for the conservation of probability. The model is defined through the rates TgigJ 
which describe the hopping process #i#2 — * I n our model the nonzero rates are the following: 



r A0 _ 

r B0 _ 7 
L 0B - 6 1 

r BA - ^ 



1 A0 - fl 2 

1 B0 - b 2 

r^-B _ _ 

1 - c 2 



The boundary terms hi and /ix, in Q act only on the first and L-th site respectively and choosing the basis 

\A) = I I . \D) 



(6) 



they have the form 
# + 7f 



-7i -«f A B + it 

where a A ( B \ and 7 + ' ' are the rates for the following processes: 
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and the indexes Z and r indicate the left and right boundary respectively. 

It is known that in the finite periodic system there is a family of steady states which are product measures 
We call P v a product measure if it has the form 

p _ p(X)p(2) p(L) 



(7) 



(8) 
(9) 
(10) 



(11) 



where P^ is a probability measure on {A, B, 0}. This means that the probability vector has the following direct 
product form: 




\P) 



where is the probability of finding an "X" particle on site k. 

Straightforward calculations show that for having a stationary product measure with uniform densities 



(12) 




(13) 
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the rates have to satisfy the following condition 

ai - a 2 - bi + b 2 + ci - c 2 = 0. (14) 

This is already a sufficient condition in case of an infinite chain or periodic boundary conditions. For open boundaries 
there are some additional restrictions for the boundary rates, for which we find: 

(ai - a 2 )p A (l p A ) (h b 2 )p A p B = a A (l-p A - p B ) (f3 A + ^)p A + 7 + p B (15) 

= [tf + 1 -)p A a A (l -p A - p B ) - j+p B (16) 

(h b 2 )p B {\ p B ) {a x a 2 )p A p B =a B (l-p A - p B ) (f3 B + 7 +)p B + jfp A (17) 

= (P B + j+)p B a B (1 - / - p B ) - j-p A . (18) 

The physical meaning of these equations is that the steady state currents of the conserved densities have to be fitted 
at the boundaries. In the following we will assume that condition l|14|) is fulfilled. For boundary rates as given above 
we say that the system is in contact with a reservoir of densities p A ' B . 

III. INVARIANT SHOCK MEASURE 

Krebs et. al. in Q pointed out that in the ASEP with open boundaries there is an invariant shock measure which 
is a product measure with a step-like density profile; p;( r ) on the left (right) of the shock. This means that if \P k ) is 
a state having the shock between the sites k and k + 1 then these states for k = 0, 1, 2, . . . , L generate a subspace of 
the vector-space of states which is invariant under time evolution and thus the many-particle problem is reduced to 
a one-particle one. This property holds if the densities pi iT - satisfy the condition 

v = PAi_jnl (19) 

q pi(l-p r ) 

where p and q are the hopping rates of the particles to the right resp. left. The resulting one-particle dynamics have 
a natural interpretation as a simple random walk of the shock position. 

In the following we search for the conditions under which there exists in our three-state model an invariant shock 
measure which is a product measure with a jump in the local particle density. To this end let \Pk) be defined as 

(20) 

The probability vectors \Pk) define an L + 1-dimensional subspace of the vector space on which H acts. Closure of 
this family of shock measures under the time evolution generated by H is equivalent to requiring 

H\P k ) = -d r \P k+1 ) - di\P k ^) + {d r + di)\P k ) for 1 < k < L - 1, (21) 
H\P ) = -d r \P x ) +d r \P ) and (22) 
H\P L ) = -di\P L -i) + di\P L ). (23) 

This is easy to see because of having only nearest neighbour interactions double hoppings of the shock position (e.g. 
\Pk) — * l-Pfc+2)) cannot occur. If this condition can be met by an appropriate choice of model parameters the quantities 
di <r and di^ r are non-negative real numbers representing the hopping rates of the shock in the bulk (di r ) and at the 
boundaries (d;,r)- 

To investigate whether I|21l23ll can be satisfied it is better to rewrite the Hamiltonian in the form 

L-l 

H = hi + (ax - a 2 )nf + (h - b 2 )n B + h kM+i + h L - (ai - a 2 )n A - (bi - b 2 )nf, (24) 

k=l 

where 

h k , k+1 = h ktk+1 - (ai - a 2 )(n A - n A +1 ) - (h - b 2 )(n B - nf +1 ). (25) 
Here and in l|24|) n A and n B are the particle number operators for the A and B particles on site k. 
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As a result hk,k+i\Pl) = for k ^ I. Therefore one has 

H\P \ = /^- fc +il Pfe ) + + hL + " a ^ n i ~ u l) + ( b i ~ fc 2)K - nf)) \P k ) if 1 < k < L- 1 
' k) \(h 1 + h L + (a 1 -a 2 ){nt-ni) + (b 1 -b 2 )(n? -nf))\P k ) iffc = 0,i 1 j 

At the boundaries one has to satisfy (|15f) with the slight modification that on the left (right) boundary we write pf' B 
(p AB ) instead of p A ' B : 

af {I -pf- P f) - {(1? + 7f)pt + !+p? = (ai - a 2 )p A (1 - pf) - (h b 2 )p A pf 

(f3 A + 7,7)^ - 7 +pf ~ a?(l - ^ - ) = (oi - a 2 )^(l - p^) - (for - Mpfpf , (27) 

and the same for the B particles: 

of (1 - Pf - Pf ) - (A B + 7, + )pf + Tfpf = (6i - h)p B (1 - pf ) - (ai - a 2 )pf pf 

(/3f + 7 +)pf - 7r"Pr - «f (1 - Pr " pf ) = (&i - Mpf (1 - pf ) - (oi - a 2 )p A p B . (28) 

This fixes a manifold of boundary parameters required to satisfy (|21I23|) . The physical interpretation is a connection 
to boundary reservoirs with different left and right densities respectively. Requiring I|27I28[) for the boundary operators 
h\ and Hl simplifies (|26[) and leads to 

H\P k ) = h k . k +i\Pk) + ((01 - a 2 )(p A - pf) + {h - 63) (pf - pf )) \P k ) for 1 < A; < L - 1. (29) 
Thus satisfying (|21|) is equivalent to requiring 

h k ,k+i\Pk) = -d r \P k +i) - di\P k -i) 

+ (dr + d l -(a 1 -a 2 )(p A ~p A )~{b 1 -b 2 )(p B -p B ))\P k ), forl<fc<L-l (30) 

Straightforward but lengthy calculation shows that hk,k+i\Pk) has this form for 1 < k < L — 1 if and only if 

ai = bi =: p, a 2 = b 2 =: q, (31) 
P _ Pr(l - Pi) 



q pi(l-p r ) 



(32) 



% = ^='-r, (33) 

Pf Pr 



where pi( r ) = pf^ + pf^ is the total density of particles on the left (right) of the shock. We note here that l|31|) 
together with (|14|l gives also c\ = c 2 . This means that the A and B particles behave the same way: the hopping rates 
to the left and to the right are p and q for both and also the exchange rates c\ and c 2 are the same in both directions. 
Equations 122I23|I give extra conditions for the boundary rates, namely: 



a t =ra l 



B and a A = ra B (34) 



This together with (|27l I28fl give conditions also for the boundary rates flf' B and 7j + r . 

The rates d r and di which describe the diffusion of the shock position to the right and to the left are 



p r l-pr , pi I -pi ,„_^ 
d r = — q = p, di = —p = - q. (35) 

Pi I— pi Pr 1 - Pr 



In addition we find the boundary rates 



d r = SI - p r (p - g) = (p - g)(l - p r ) + r f\^ y (36) 
Pi (r + l)(l-pj) 

dl = ^ + Pl ( P -q) = -(? - g)(l - Pi) + , ^tt^ y (37) 

p r (r+l)(l-p r ) 



where a l(r) = a A [r) + a B r) . 
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Using (|32|) the bulk shock hopping rates (|35|l become 

, , s Pl(x)0- - Pl{r)) i , 

di(r) = (p - q) (38) 

Pr ~ Pi 

as in the ASEP 0, . Since by identifying A and B particles one arrives at the one-species ASEP, it is not surprising 
that ll.')2[l and l|38[) are in full agreement with the corresponding formulas of the ASEP. The novelty is that a.) 
the product measure with step-like density profiles satisfying l)33|) remains an invariant shock measure even if one 
distinguishes between A and B particles b.) there is no other possibility for having such an invariant shock measure 
(up to relabelling the states A, B and 0). Notice that this solution describes a single shock whereas generically a 
system with two conservation is expected to have two stationary shocks 



IV. HYDRODYNAMIC LIMIT 
A. Derivation of the partial differential equation 

In the previous sections the microscopic dynamics of the two-species ASEP was studied. It was found that invariant 
product measures with step- like density profiles are possible if criterion i|31|) for the rates and criteria Q32p . I|33|) for 
the densities are fulfilled. It is also of interest to drop the conditions (|32fl . (|33|l and see how this model behaves on 
the macroscopic scale, i.e., when the space and time is rescaled as i mac = t m - lc a, x — ka and the lattice spacing a — > 
(Eulerian scaling). In this subsection we restrict the discussion to the case of infinite chains (or torus geometry). 
Finite systems with open boundaries are considered in the last subsection. 

Performing the Euler scaling one can easily construct the naive hydrodynamic equations of a model having two 
conserved quantities (v and u) and stationary product measure. These describe the macroscopic time-evolution of 
the conserved quantities: 

d t v + d x j v (v,u) = (39) 
d t u + d x j u {v,u) = 0, (40) 

where j v (v,u) and j u {v,u) are the currents of v and u assuming homogeneous product measure. We have to note 
here that it is far from trivial that these set of conservation laws describe correctly the macroscopic dynamics of the 
model. However, a mathematically rigorous proof is available that in all systems having nearest neighbor dynamics 
with two conserved densities (v and u) and stationary product measure, the time evolution of the densities under 
Eulerian scaling is described by th e eq uations H39I40J) until the occurrence of the first shock [8j. For closely related 
results on two-species systems, see |l3j . for an apparent failure of the hydrodynamic description, see [l4| and earlier 
work on condensation in two-species systems |lCt 1 151. Il6j . 

For our model it is useful to introduce the new set of conserved quantities (v,u) instead of (p A ,p B ): 

v = 1 — p A — p B (the density of vacancies) (41) 

u = p A ~p B . (42) 

The associated currents are: 

j v = -(p-q)v(l-v) (43) 

ju = {p- q)vu. (44) 

Plugging these into (|39I4U[1 one arrives to the hydrodynamic equations of the model: 

d t v - (p - q)d x (v(l - v)) = (45) 
d t u + {p - q)d x {vu) =0. (46) 

In the special case when p = q the currents are zero and the solution of (|39I4U|I becomes trivial. In this case the 
behaviour of the microscopic model is diffusive so one has to perform diffusive scaling (i mac = t m i c a 2 , x = ka, a — > 0) 
to see the nontrivial macroscopic behaviour. 

It is to be noted that the generic set of PDEs describing systems with three local states and satisfying l|14|l is the 
so-called Leroux's system 

d t a + d x {ar) = (47) 
d t T + d x (a + r 2 ) = 0. (48) 
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It means that there are always such linear combinations (a, r) of the conserved quantities (shifted by an irrelevant 
constant) which satisfy (|47I48|) 8J. However, in our case, when a\ = b\ and a 2 = b 2 , and r turn out to be linearly 
dependent and equations 147(1 and 1(48(1 are equivalent (and correspond to 1(45(1 ) . Therefore our model has to be 
investigated separately. 

Note that the equation 145(1 for v is decoupled from u and takes the form of the well-known Burgers equation, which 
can be solved exactly 0|- Introducing 

cb=(p A -p B )/(p A + p B ) =u/(l-v) (49) 

one gets the linear wave equation for <fi: 

d t <t>+{p-q)vd x <t> = (50) 
meaning that </> is constant along the curves x(t) satisfying the ordinary differential equation: 

^=v(x,t), (51) 

where v(x, t) is the solution of ((45(1 . The physical meaning of x(t) is that this is the expected path of a single tagged 
particle (either A or B). 



B. Development of an initial sharp interface 



Suppose that at t = there is a sharp interface separating domains characterised by u;( r ) and uim on the left 
(right) respectively. We assume that vi > v r in order to guarantee the stability of the interface No further 

conditions on the values of u, v are imposed. 

The time and space dependence of v is known since it is described by the inviscid Burgers equation 1(45(1 : the 
step-like front (having on the left (right) of it) is travelling with velocity 

V 1 = (p-q){l-v l -v r ). (52) 

The x(t) curves defined by equation 1(51(1 are then broken lines as in figure^ One can see that (j) changes while crossing 
the bold dashed line indicating that there is another shock in the system travelling with larger velocity (without loss 
of generality we assume that p > q) 

V 2 = {p-q)v l >V 1 . (53) 

Generally a single sharp interface at t ~ forms two shocks (or rarefaction waves) at t > in systems having two 
conservation laws |6( . For the complete analysis one has to study the eigenvalues and eigenvectors of the Jacobian 

d{juJv) -(p-<l)( V no„\) (54) 



d(u,v) * HJ V 2v ~ 1 

which give the characteristic velocities of the different types of density fluctuations . In this specific system these 
are: 

\ 1 = (p-q)(2v-l) ex^^j) (55) 

A 2 = ip - q)v e 2 = y q ^ (56) 

The two corresponding types of shocks have either constant cj) across the discontinuity (type 1, occurring along the 
solid bold line in Fig. 1) or constant v (type 2, occurring along the broken line in Fig. 1). A shock (microscopically 
sharp interface) satisfying either of these relations initially does not split, whereas a generic shock splits into two 
shocks of type 1 and type 2 respectively. We remark that condition 1(33(1 translates into tfii — <\> 2 for the quantity </>. 
Hence the shock studied in the previous section is of type 1. Notice also that in the second mode u is changed but v 
is not. Walking along these lines in the (u, v) space Ai remains unchanged. This means that if there is a discontinuity 
in u (with constant v) then the fluctuations in u are not driven towards the "shock" but they travel with the same 
mean velocity resulting in an unstable shock on the microscopic scale. However, since in the original lattice model 



x 



X 



FIG. 1: The x(t) curves of equation (^J for the case when v(x, t) is a shock-solution of the Burgers equation ■ The position 
of this shock (starting from xo) is marked by the thick solid line, u also changes on the dashed line indicating another shock 
in the system, but this is not microscopically sharp. 

the width of the step is expected to scale with y/i due to the diffusive dynamics of A and B particles the interface 
remains sharp under Eulerian scaling. 

The physical significance of two types of shocks can be easily demonstrated considering the following example: at 
t = the system is partially occupied with only A particles on the left of xq and fully occupied by B particles on 
the right. This setting models a situation where particles (A) fall down (e.g. due to gravity) and enter a medium 
(B) where they have no weight. Hence they penetrate diffusively (with vanishing macroscopic flux), but pile up 
ballistically due to the constant incoming macroscopic flux. The result are two interfaces emerging (see figure [2}. The 
second (lower) interface has velocity since v\ = V2 = 0, see (|53|l . Although the particles diffuse and start mixing on 
the lattice scale at the A — B interface, this remains invisible on the Euler scale. 

The other interface (type 1) separating domains with different densities of A particles corresponds to the usual 
shocks which are known from the study of the Burgers equation. Here both u and v are changed (but (j) remains 
constant). 

On figure |3 we show a more general example for the time development of the system starting from an initial state 
having a step in both densities at the origin. The second mode also can develop rarefaction waves if the density on 
the left is lower than on the right in the initial state. We note that the shock coming from the second mode always 
stays to the right of the shock/rarefaction wave coming from the first mode since it is faster 11531 . This is by the 
relation v > 2v — 1 (rarefaction waves) and V1.2 > vi + V2 — 1 (shocks). 

C. Steady state selection in the open system 

In a finite system with open boundaries one is left with the question of steady state selection. In the infinite 
system all the states with constant density profiles are stationary, however in the finite system with boundaries acting 
as particle reservoirs the bulk densities depend non-trivially on the boundary densities. Setting d t = in |JS3J||@DJ| 
one arrives at the solution that u and v are constant, so in general the densities cannot fit both boundaries. This 
discrepancy is resolved by the appearance of shocks in the original lattice model leading to discontinuities at either 
(or both) boundary in the hydrodynamic limit. 

In case of one conserved density the density-current relation already determines the phase diagram in terms of the 
boundary densities |l8l [l9l | but in case of more conservation laws the question turns out to be much more intricate 
and no general rule is known to apply [7j. However, in our model it is possible to determine the resulting steady state 
for given boundary densities. 

Since the dynamics of v is i ndep endent of u and follows the usual ASEP dynamics the profile can be deduced from 
the ASEP phase diagram [2(J, |2j, [22 • Our task is then to determine the bulk value of u in terms of the boundary 



t = 



t > 





FIG. 2: The time development from a special initial condition: A particles fall down (e.g. due to gravity) and enter a medium 
where they have no weight, hence they penetrate diffusively but pile up ballistically resulting in two interfaces. The width of the 
domain wall between the A (black) and B (gray) particles scales with \ft while the other domain wall remains microscopically 
sharp (however, on a diffusive scale one could see the fluctuations of the position of this sharp interface). 
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(v 1 +v 2 -l)t 



Vot 



FIG. 3: Development of the densities v and u under Eulerian scaling starting from an initial state with step- like density profiles: 
Vi,ui on the left and V2,U2 on the right of the origin. The first interface (on the left) is a shock coming from the Burgers- 
equation and corresponds to the thick solid line on figure The second one (on the right) is not microscopically sharp. On 
figure it is the dashed line. 



densities u lclt , u light , v\ clt , v light . 

Since all the fluctuations in u travel with velocity v > one would think that itbuik is only determined by the left 
boundary if v ^ 0. However, this is not the case because a discontinuity in v can be localised at this boundary which 
also induce a step in u. The previously introduced </>, however, does not change at this discontinuity which suggests 
that the bulk value of </> is given by the left boundary (0buik = 0ioft)- In order to verify this heuristic reasoning we 
make use of a technique to determine the steady state which is widely used for systems with a single conservation 
law. One introduces a viscosity term in I|39I40|I proportional to v which contain a second derivative of the densities 
leading to second order equations for the stationary profiles. Now the profiles can fit both boundaries for any value 
of v and after performing limit v — > one can get the selected density. 

In case of one single conservation law this technique is quite robust in the sense that the resulting bulk density is 
essentially independent of the specific choice of the viscosity term (as long as it contains a second order derivative) . 
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V left I/ 2 




1/2 



V 



right 



FIG. 4: The phase diagram of the open system. Since the vacancies follow the ASEP dynamics the boundary values of v 
determine the phase according to the ASEP phase diagram. In addition equation (1571 determines the stationary value of u in 
all phases except on the dashed line where one has symmetric diffusion of A and B resulting in linear profiles for it connecting 

Mleft/(1 - Wlcft) and Wright ■ 



For more conservation laws this method is mathematically not well-studied and one has to be more careful with the 
choice of viscosity term. For systems with a stationary product measure there is a natural choice for these terms 
introduced in Q. Namely, when one derives the hydrodynamic equations from the lattice continuity equations one 
uses a Taylor expansion of the current (as a function of the densities) in the lattice constant. Keeping the second 
order terms leads to a unique viscosity term which is proportional to the lattice constant a. 

In appendix^ we show that using this technique for our model the above conjecture, that </>buik = </>i e ft, is confirmed 
if Wright 7^ 0, which implies 

1 - v 

U = Uleft- if Wright ^ 0- (57) 

1 - Uleft 

Summarizing the results (see also figure 2J) we find 

{V ie ft if Uleft > 1/2 and U righ t > 1 - Uleft 
Wright if "right < 1/2 and v light < 1 - uieft , (58) 
1/2 if Uieft < 1/2 and u rlght > 1/2 

^ f Uleft if Uleft > 1/2 and Urfght > 1 - Uleft 

"bulk = Mleft , _ bUlk = I "left ^-vlltt if < Wri S ht - l l 2 and Wri S ht < 1 ~ Vlelt ■ ( 59 ) 

' 1 [ Mieft i-^Lt if vieft - 1//2 and Viisht - 

On the coexistence line (0 < u r i g ht = 1 — ui e ft < 1/2) the stationary w-profile is linear as known from the ASEP 
pfjl l2ll| . This induces also a linear profile for u according to (|57|l . 

If Ujight = then there are no vacancies in the system leading to symmetric diffusion of A and B particles. This 
results in linear profiles for u connecting ui c ft/(l — ui e ft), which is the value of u after the discontinuity located at the 
left boundary, and w r i g ht as expected from the symmetric exclusion process (SEP) |23l | and confirmed in appendix lAl 
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V. CONCLUSION 

We have introduced an asymmetric exclusion process with two conserved species of particles. Its study is part 
of the program to understand the emergence of macroscopic behaviour from microscopic models of nonequilibrium 
diffusive interacting particles and to elucidate the significance of boundary conditions in driven diffusive systems. 
Under heuristic Eulerian scaling we obtain a system of PDE's which can be analysed analytically. As one would have 
expected from the nonlinear current-density relation the system may produce shock discontinuities on the macroscopic 
scale. There are two types of shocks (as expected in systems with two conservation laws) for which we calculate the 
respective shock velocities and intermediate densities in terms of the asymptotic limiting densities. 

The model we have introduced here allows for a rigorous analysis of the microscopic structure of a shock of type 1 
(with a discontinuity in the vacancy density), provided that some special conditions on the hopping parameters and 
shock densities are satisfied. For this shock one finds the same behaviour as for the ASEP, suggesting that also for 
nongeneric values of the limiting densities of the shock the behaviour would be similar. Shocks of type 2 (constant 
vacancy density) are microscopically unstable, they correspond to diffusive spreading as in the symmetric exclusion 
process. 

Part of our investigation is devoted to the problem of steady state selection in the open system [24] . Unlike for 
driven diffusive system with a single species of particles there is no developed theory of boundary-induced phase 
transitions for systems with two or more conservation laws. Using the known phase diagram of the ASEP and the 
microscopic properties of the shock in our particular model we have derived the stationary phase diagram. So our 
model, even though degenerate, may serve as a testing ground for any general theory. With a view on some of our 
earlier results on driven diffusive two-species systems with open boundaries we note that we found that the heuristic 
approach to steady-state selection of Ref. jfj reproduces the independently derived phase diagram. 

An intriguing observation is that the u-density in the phase diagram depends continuously on the left w-density 
throughout the maximal current phase. In order to explain the significance of this observation we note that the 
phase diagram of single-species systems can be predicted from the current-density relation in terms of boundary 
densities [L&L Il9| , which in general are unknown functions of the boundary rates of the model. On the other hand 
any numerical study of a given model yields the phase diagram in terms of the boundary rates, not the boundary 
densities. Therefore a comparison between theoretical predictions and numerical observations requires knowledge of 
the nonuniversal relationship between boundary densities and boundary rates. In the ASEP and generally in single- 
species systems the order parameter (the bulk density) in the maximal current phase does not depend on the left 
boundary density. Hence it is not possible to measure the relationship between the left boundary density and the 
boundary rates of the model. This in turn makes it impossible to predict the phase diagram in terms of the boundary 
rates. In the present two-species model, however, the left boundary density can be measured in the maximal current 
phase as a function of the system parameters through measuring u. 

This leads to an unexpected offspin of our investigation: by regarding the two-species model as an exclusion 
process with two kinds of vacancies (by identifying vacancies with particles and the two kinds of particles with 
vacancies). Generalizing other single-species systems in a manner similar to what is done here, viz., leaving their 
dynamics unchanged by just introducing tagged diffusive vacancies, may provide a means of measuring those postulated 
boundary densities in terms of the model parameters even in parts of the phase diagram where the order parameter 
itself does not depend on the boundary densities. 

Finally we remark that in principle the processes A — > B or B — > A could also be excluded (i.e. liiyf* = 0) 
from our process. We only keep them for the sake of generality. If one wants to choose the boundary conditions to 
mimic certain boundary densities these rates are in general non-zero. However, in our model these processes don't 
play an important role unlike in the so-called bridge model where they are responsible for spontaneous symmetry 
breaking |25| . The precise mechanism for this phenomenon and a quantitative description, however, still remain a 
major challenge in the study of two-species systems with open boundaries. 
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APPENDIX A: SOLUTION OF THE HYDRODYNAMIC EQUATIONS WITH FINITE VISCOSITY 

The viscous hydro-dynamic equations for the steady state coming from the method described in Ref. @ (see Sec. 
IIV Cp are the following: 

(p-q) v >(2v-l) = a^-v", (Al) 
(p - q)(u'v + v'u) = cJ^-{u"v - v"u) + ac ((1 - v)u" + uv") . (A2) 
Substituting u = (1 — v)(j> in l|A2() and using IjAlfl we arrive to the following equation for <f>: 

v h _ V U> = a (P + l) f v n _ w u» _ 2vv' ( f ) ') + — ((1 - w) V - 2(1 - vU'v') . (A3) 
2{p - q) P~q 

One can immediately see here that <j) — const, is a solution if it hts both boundaries, i.e., if 4>\ c n — bright- When this 
does not hold (this is the case in general) then one has 

— In \<b'\ = 1 — r r . (A4) 

Integrating from the right boundary we get 

Assuming that the sign of <fi' does not change in the system we arrive to 

For strictly positive v in the bulk of the system one has <fi'(y) — > as a — > for any fixed y < 1 thus the limiting <p 
profile is flat everywhere apart from the right boundary, which implies 0tmik = <^icft- 

There is no positive lower bound for v if and only if u r i g ht = (dashed line in figure 0J- In this case v(x) — > 
everywhere as a — > apart from the left boundary (provided Ui e ft ^ 1) which implies 4>"(x) — > for any fixed x > 
according to (|A3() . We still have to show that <fi does not have a discontinuity in the a — > limit at the left boundary. 
For this we evaluate <fi'(0) and show that it does not diverge. Similarly to l|A6|l one has 

0'(O) =0U f-^) 2 exp (-^ /" u -dx j , (A7) 

where the upper limit of the integration < y < 1 is arbitrary because of the fast decay of the integrand. Since the 
integrand is non-negative the exponential cannot diverge which gives 

lim0'(O)<oo. (A8) 

This means that the <p profile is the linear function connecting the two boundary values. The u profile, however, has a 
discontinuity at the left boundary according to (|49|l and jumps from Ui e ft to ui c ft/(l — vicit)- From here it goes linearly 
to the right boundary value M r i g ht- 
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